Geodetic mass balance record with rigorous uncertainty estimates deduced from aerial photographs and lidar data – Case study from Drangajökull ice cap, NW Iceland

In this paper we describe how recent high-resolution digital elevation models (DEMs) can be used to extract glacier surface DEMs from old aerial photographs and to evaluate the uncertainty of the mass balance record derived from the DEMs. We present a case study for Drangajokull ice cap, NW Iceland. This ice cap covered an area of 144 km(2) when it was surveyed with airborne lidar in 2011. Aerial photographs spanning all or most of the ice cap are available from survey flights in 1946, 1960, 1975, 1985, 1994 and 2005. All ground control points used to constrain the orientation of the aerial photographs were obtained from the high-resolution lidar DEM. The lidar DEM was also used to estimate errors of the extracted photogrammetric DEMs in ice-and snow-free areas, at nunataks and outside the glacier margin. The derived errors of each DEM were used to constrain a spherical semivariogram model, which along with the derived errors in ice-and snow-free areas were used as inputs into 1000 sequential Gaussian simulations (SGSims). The simulations were used to estimate the possible bias in the entire glaciated part of the DEM and the 95% confidence level of this bias. This results in bias correction varying in magnitude between 0.03m (in 1975) and 1.66m (in 1946) and uncertainty values between +/- 0.21m (in 2005) and +/- 1.58m (in 1946). Error estimation methods based on more simple proxies would typically yield 2-4 times larger error estimates. The aerial photographs used were acquired between late June and early October. An additional seasonal bias correction was therefore estimated using a degree-day model to obtain the volume change between the start of 2 glaciological years (1 October). This correction was largest for the 1960 DEM, corresponding to an average elevation change of -3.5m or approx. three-quarters of the volume change between the 1960 and the 1975 DEMs. The total uncertainty of the derived mass balance record is dominated by uncertainty in the volume changes caused by uncertainties of the SGSim bias correction, the seasonal bias correction and the interpolation of glacier surface where data are lacking. The record shows a glacier-wide mass balance rate of (B) over dot = -0.26 +/- 0.04m w.e.a(-1) for the entire study period (1946-2011). We observe significant decadal variability including periods of mass gain, peaking in 1985-1994 with (B) over dot = -0.27 +/- 0.11m w.e.a(-1). There is a striking difference when (B) over dot is calculated separately for the western and eastern halves of Drangajokull, with a reduction of eastern part on average similar to 3 times faster than the western part. Our study emphasizes the need for applying rigorous geostatistical methods for obtaining uncertainty estimates of geodetic mass balance, the importance of seasonal corrections of DEMs from glaciers with high mass turnover and the risk of extrapolating mass balance record from one glacier to another even over short distances.This work was carried out within SVALI funded by the Nordic Top-level Research Initiative (TRI) and is SVALI publication number 70. It was also financially supported by alpS GmbH. This work is a contribution to the Rannis grant of excellence project, ANATILS. We thank the National Land Survey of Iceland and Loftmyndir ehf. for acquisition and scanning of the aerial photographs. This study used the recent lidar mapping of the glaciers in Iceland that was funded by the Icelandic Research Fund, the Landsvirkjun Research Fund, the Icelandic Road Administration, the Reykjavik Energy Environmental and Energy Research Fund, the Klima- og Luftgruppen (KoL) research fund of the Nordic Council of Ministers, the Vatnajokull National Park, the organization Friends of Vatnajokull, the National Land Survey of Iceland and the Icelandic Meteorological Office.Peer Reviewe

The geodetic mass balance of Eyjafjallajökull ice cap for 1945–2014: processing guidelines and relation to climate

Publisher's version (útgefin grein)Mass-balance measurements of Icelandic glaciers are sparse through the 20th century. However, the large archive of stereo images available allows estimates of glacier-wide mass balance in decadal time steps since 1945. Combined with climate records, they provide further insight into glacier-climate relationship. This study presents a workflow to process aerial photographs (1945-1995), spy satellite imagery (1977-1980) and modern satellite stereo images (since 2000) using photogrammetric techniques and robust statistics in a highly automated, open-source pipeline to retrieve seasonally corrected, decadal glacier-wide geodetic mass balances. In our test area, Eyjafjallajökull (S-Iceland, ~70 km2), we obtain a mass balance of <![CDATA[$, with a maximum and minimum of and , respectively, attributed to climatic forcing, and , mostly caused by the April 2010 eruption. The reference-surface mass balances correlate with summer temperature and winter precipitation, and linear regression accounts for 80% of the mass-balance variability, yielding a static sensitivity of mass balance to summer temperature and winter precipitation of-2.1 ± 0.4 m w.e.a-1K-1 and 0.5 ± 0.3 m w.e.a-1 (10%)-1, respectively. This study serves as a template that can be used to estimate the mass-balance changes and glaciers' response to climate.This study was funded by the University of Iceland (UI) Research Fund. Collaboration and travels between IES and LEGOS were funded by the Jules Verne research fund. We thank David Shean and two anonymous reviewers for their valuable comments, which greatly improved the manuscript. We thank Carsten Kristinsson at LMÍ for scanning the aerial photographs, Oleg Alexandrov for his helpful tips and advice on ASP, Luc Girod for his help in the MicMac forum and Deirdre Clark and Ken Moxham for the Englishlanguage editing of the manuscript. Pléiades images were acquired at research price thanks to the CNES ISIS program (http://www.isis-cnes.fr). This study uses the lidar mapping of the glaciers in Iceland, funded by the Icelandic Research Fund, the Landsvirkjun research fund, the Icelandic Road Administration, the Reykjavík Energy Environmental and Energy Research Fund, the Klima-og Luftgruppen research fund of the Nordic Council of Ministers, the Vatnajökull National Park, the organization Friends of Vatnajökull, LMÍ, IMO and the UI research fund. This study uses the GLIMS database of the outlines of Icelandic glaciers. E.B. acknowledges support from the French Space Agency (CNES) through the TOSCA program.Peer Reviewe

Non-surface mass balance of glaciers in Iceland

Publisher's version (útgefin grein)Non-surface mass balance is non-negligible for glaciers in Iceland. Several Icelandic glaciers are in the neo-volcanic zone where a combination of geothermal activity, volcanic eruptions and geothermal heat flux much higher than the global average lead to basal melting close to 150 mm w.e. a−1 for the Mýrdalsjökull ice cap and 75 mm w.e. a−1 for the largest ice cap, Vatnajökull. Energy dissipation in the flow of water and ice is also rather large for the high-precipitation, temperate glaciers of Iceland resulting in internal and basal melting of 20–150 mm w.e. a−1. The total non-surface melting of glaciers in Iceland in 1995–2019 was 45–375 mm w.e. a−1 on average for the main ice caps, and was largest for Mýrdalsjökull, the south side of Vatnajökull and Eyjafjallajökull. Geothermal melting, volcanic eruptions and the energy dissipation in the flow of water and ice, as well as calving, all contribute, and thus these components should be considered in mass-balance studies. For comparison, the average mass balance of glaciers in Iceland since 1995 is −500 to −1500 mm w.e. a−1. The non-surface mass balance corresponds to a total runoff contribution of 2.1 km3 a−1 of water from Iceland.Financial support for lidar mapping of glaciers in Iceland in 2008–2012 was provided by the Icelandic Research Fund (163391-052), the Landsvirkjun (National Power Company of Iceland) Research Fund, the Icelandic Road Administration, the Reykjavík Energy Environmental and Energy Research Fund, the National Land Survey of Iceland, the Klima- og Luftgruppen (KoL) research fund of the Nordic Council of Ministers, and the Vatnajökull National Park. The acquisition of the Hofsjökull 2013 DEM was funded by AlpS GmbH and the University of Innsbruck. The acquisition of the Langjökull 2013 DEM was funded by NERC grant IG 13/12 and the DEM was provided by Ian Willis at the Scott Polar Research Institute. The work on estimating geothermal and volcanic power is based on funding from many sources, including the Research Fund of the University of Iceland, ISAVIA (the Icelandic Aviation Service), the Icelandic Road Administration and Landsvirkjun; logistical support has been provided by the Iceland Glaciological Society.Peer Reviewe

Hekla Volcano, Iceland, in the 20th Century: Lava Volumes, Production Rates, and Effusion Rates

Publisher's version (útgefin grein)Lava flow thicknesses, volumes, and effusion rates provide essential information for understanding the behavior of eruptions and their associated deformation signals. Preeruption and posteruption elevation models were generated from historical stereo photographs to produce the lava flow thickness maps for the last five eruptions at Hekla volcano, Iceland. These results provide precise estimation of lava bulk volumes: V1947–1948 = 0.742 ± 0.138 km3, V1970 = 0.205 ± 0.012 km3, V1980–1981 = 0.169 ± 0.016 km3, V1991 = 0.241 ± 0.019 km3, and V2000 = 0.095 ± 0.005 km3 and reveal variable production rate through the 20th century. These new volumes improve the linear correlation between erupted volume and coeruption tilt change, indicating that tilt may be used to determine eruption volume. During eruptions the active vents migrate 325–480 m downhill, suggesting rough excess pressures of 8–12 MPa and that the gradient of this excess pressure increases from 0.4 to 11 Pa s−1 during the 20th century. We suggest that this is related to increased resistance along the eruptive conduit.Icelandic Research Fund. Grant Number: 152266‐052Peer Reviewe

Lava field evolution and emplacement dynamics of the 2014–2015 basaltic fissure eruption at Holuhraun, Iceland

The 6-month long eruption at Holuhraun (August 2014–February 2015) in the Bárðarbunga-Veiðivötn volcanic system was the largest effusive eruption in Iceland since the 1783–1784 CE Laki eruption. The lava flow field covered ~84 km2 and has an estimated bulk (i.e., including vesicles) volume of ~1.44 km3. The eruption had an average discharge rate of ~90 m3/s making it the longest effusive eruption in modern times to sustain such high average flux. The first phase of the eruption (August 31, 2014 to mid-October 2014) had a discharge rate of ~350 to 100 m3/s and was typified by lava transport via open channels and the formation of four lava flows, no. 1–4,which were emplaced side by side. The eruption began on a 1.8 km long fissure, feeding partly incandescent sheets of slabby pāhoehoe up to 500 m wide. By the following day the lava transport got confined to open channels and the dominant lava morphology changed to rubbly pāhoehoe and ‘a’ā. The latter became the dominating morphology of lava flows no. 1–8. The second phase of the eruption (Mid-October to end November) had a discharge of ~100–50 m3/s. During this time the lava transport system changed, via the formation of a b1 km2 lava pond ~1 km east of the vent. The pond most likely formed in a topographical low created by a the pre-existing Holuhraun and the newHoluhraun lava flow fields. This pond became themain point of lava distribution, controlling the emplacement of subsequent flows (i.e. no. 5–8). Towards the end of this phase inflation plateaus developed in lava flowno. 1. These inflation plateaus were the surface manifestation of a growing lava tube system, which formed as lava ponded in the open lava channels creating sufficient lavastatic pressure in the fluid lava to lift the roof of the lava channels. This allowed new lava into the previously active lava channel lifting the channel roof via inflation. The final (third) phase, lasting from December to end-February 2015 had a mean discharge rate of ~50 m3/s. In this phase the lava transport was mainly confined to lava tubes within lava flows no. 1–2, which fed breakouts that resurfaced N19 km2 of the flow field. The primary lava morphology from this phase was spiny pāhoehoe, which superimposed on the ‘a’ā lava flows no. 1–3 and extended the entire length of the flow field (i.e. 17 km). Thismade the 2014–2015 Holuhraun a paired flow field,where both lava morphologies had similar length. We suggest that the similar length is a consequence of the pāhoehoe is fed from the tube systemutilizing the existing ‘a’ā lava channels, and thereby are controlled by the initial length of the ‘a’ā flows.The work was financed with crisis response funding from the Icelandic Government along with European Community's Seventh Framework Programme Grant No. 308377 (Project FUTUREVOLC) and along with the Icelandic Research fund, Rannis, Grant of Excellence No. 152266-052 (Project EMMIRS). Furthermore, Vinur Vatnajökuls are thanked for support.Peer Reviewe

Integration of SAR Data Into Monitoring of the 2014-2015 Holuhraun Eruption, Iceland: Contribution of the Icelandic Volcanoes Supersite and the FutureVolc Projects

We report how data from satellite and aerial synthetic aperture radar (SAR) observations were integrated into monitoring of the 2014–2015 Holuhraun eruption in the Bárðarbunga volcanic system, the largest effusive eruption in Iceland since the 1783–84 Laki eruption. A lava field formed in one of the most remote areas in Iceland, after the propagation of a ∼50 km-long dyke beneath the Vatnajökull ice cap, where the Bárðarbunga caldera is located. Due to the 6 month duration of the eruption, mainly in wintertime, daily monitoring was particularly challenging. During the eruption, the European volcanological project FutureVolc was ongoing, allowing collaboration of many European experts on volcano monitoring activities. Icelandic volcanoes are also a permanent Supersite within the Geohazard Supersites and Natural Laboratories (GSNL) initiative, with support from the Committee on Earth Observation Satellite (CEOS) in the form of a large collection of SAR images. SAR data were acquired by Cosmo-SkyMed (CSK) and TerraSAR-X (TSX) satellites and complemented by aerial SAR images. The large set of SAR satellite data significantly contributed to the daily monitoring during the unrest at Bárðarbunga caldera, the Holuhraun eruption and the year following the eruption. Detection of surface changes using both SAR amplitude and phase information was conducted throughout the whole duration of the volcano-tectonic event, and in the following months, to quantify and track the evolution of volcanic processes at Holuhraun and geothermal activity at Bárðarbunga volcano. Combination of SAR data with other data sets, e.g., satellite optical images and geodetic Global Positioning System (GPS) measurements, was essential for the evaluation of the volcanic hazard in the whole area. International collaboration within the FutureVolc project formed the basis for successful analyses and interpretation of the large SAR data set. Information was provided at Scientific Advisory Board meetings of the Icelandic Civil Protection and used in decision-making, as well as for supporting field-deployment and air-based surveys